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Accurate numerical solution of the five-body Schrodinger equation is effected via variational Monte 
Carlo. The spectrum is assumed to exhibit a narrow resonance with strangeness S = +1. A fully 
antisymmetrized and pair-correlated five-quark wave function is obtained for the assumed non- 
relativistic Hamiltonian which has spin, isospin, and color dependent pair interactions and many- 
body confining terms which are fixed by the non-exotic spectra. Gauge field dynamics are modeled 
via flux tube exchange factors. The energy determined for the ground states with J'^ = + ) is 
2.22 GeV (2.50 GeV). A lower energy negative parity state is consistent with recent lattice results. 
The short-range structure of the state is analyzed via its diquark content. 



A system of interacting, non-relativistic constituent 
quarks is the most simple, realistic model of hadronic 
systems. Solving the many-body Schrodinger equation 
to determine wave functions within this simple model is 
still a formidable task owing to the strong flavor, spin, 
and color dependence of the quark-quark interaction and 
traditionally requires some array of approximate meth- 
ods to solve it. The controversial status of the recent 
evidence of a flavor exotic five-quark state warrants a 
careful treatment of this strongly interacting system. In 
this letter, the exercise of determining the wave function 
of five interacting constituent quarks in the fiavor-exotic 
multiquark hadronic sector is solved using the variational 
Monte Carlo (VMC) technique. This technique is known 
to yield upper bounds on the ground state energy accu- 
rate to the level of a few percent in light nuclei with the 
number of nucleons A < 6 . 

Recent experimental evidence 
Ell IT3. has revived interest in the multiquark flavor- 
exotic sector of the hadronic spectrum. Various model 
calculations have been reported [H, [H, EE E3, E3, E 
which study the existence of a strangeness 
S = +1 resonance, dubbed 9'^, about 100 MeV above 
threshold to nucleon-kaon decay with low mass of 1540(2) 
MeV and possibly extremely narrow width of 0.9(3) MeV 
^1^ . Lattice results are available, including Refs. [l^ 

The calculation of Ref . , working in the chiral quark 
soliton model, of an extremely small width, < 15 MeV 
stimulated experimental searches with some positive sig- 
nals, although all of them had very low statistics with at 
most ~ 100 counts above a comparable background 0|. 
Preliminary high statistics photoproduction data of the 
reaction nKs from the CLAS collaboration 

sets an upper limit on the yield of the relative to the 
A* (1520) yield at 0.2%. No definitive structure in the 
ni^T-invariant mass spectrum is observed at 1540 MeV in 
this experiment. No experimental information is avail- 
able on the spin or parity of the state in any experiment 
done to date. 

We study the 6*+ in the non-relativistic flux tube quark 



model with one-gluon exchange (OGE) and one-pion ex- 
change (OPE) (between the light quarks). We take se- 
riously the possibility that 6^ is a narrow resonance 
and calculate its mass as a stable state with respect to 
strong interactions. The approach adopted here will be 
to work with a general, completely antisymmetric wave 
function and determine dynamically which flavor-spin- 
color-orbital (TSCL) structures are favored in the con- 
stituent quark model (CQM). The model Hamiltonian 
used in this work is fixed by the single hadron spectrum 
[iij and the six-quark (deuteron) properties 29]. It is 
given as a sum of kinetic energy, pair potentials for OGE 
and OPE, and confining terms: 
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Here, rrii is the constituent quark mass, — n — Yj, 
and ^/a is the flux-tube tension. The two-body opera- 
tor potentials v'^j^\rij) for OGE (OPE) interactions are 
determined in the non-relativistic reduction of the tree 
level amplitude for qq or qq scattering. The many-body 
confining interaction 

V^^iR) = V^L(R) - i ^ V^r., (2) 

i<j=i 

where R = (ri , . . . , rs) and i(R) is the total length of the 
fiux tubes, is defined with respect to a two-body contribu- 
tion (the second term) which ensures that ^■'^"^(R) 
when one of the quarks is pulled far from the others. The 
locations of the three Y junctions, denoted by (g) in the 
diagrams in Fig.(Q are determined by minimizing the 
total length of the flux tubes given the positions R of 
the quarks and topology of the tubes. Possible flux tube 
topologies for the Aqq states consistent with gauge invari- 
ance are shown in Figs.l^. They correspond to overall 
singlet states of the form 
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where ijk = 234, 432, or 324. Off-diagonal matrix ele- 
ments of operators with color dependence (like those in 
v^{rij)) of these orthonormal states are accompanied by 
a factor which models the dynamics of the gauge fields: 



(4) 



where is the range over which flux tubes can be 
exchanged (roughly the width of the flux tube ) an d rn 
is the distance between the exchanged quarks '29'|. We 
use jFT — 0, which has been determined to minimize 
the energy. The two-body confining interaction jy/arij 
of Eq.lQJ is accounted for in VMC non-perturbatively 
while the V^^ , typically ~ 15% of the total confinement 
contribution, is evaluated perturbatively. The constant 
Vb(4qq) in Eq.Q, extrapolated from the single hadron 
spectrum, will be discussed below. 

The state vector for the five quarks, including two- 
body isospin, spin, color and spatial correlations is given 

by 
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where the symmetrized product of two-body correlation 
operators acts on the uncorrelated TSCL state: 

c=R..G,B 

X [|(T5Ci)[i4,; T,, S,,c, Ml) |s; s., c, me)] 



(6) 

Here we have four light quarks, uudd (or Aq) described by 
the state |(rS'CL)[i4];T,,S'^,c,ML) with isospin, T,T^, 
spin S, Sz , color C, c, and orbital angular momentum 
L, Ml coupled to the completely antisymmetric state, 
denoted {TSCL)[n]. The flavor is assumed to be T = 
and total spin of 4g can takes the values S = 0,1, 2. We 
neglect the states with 5 = 2 since they are expected 
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FIG. 1: Topological configurations of flux tubes. In general, 
configurations of quarks do not lie in the plane. 



|i;i+> = 

12; 1+) = 
l3;(0,i)+>= 
|4; (0,1)+) = 



TABLE I: Inner products on S4 of spin, isospin, color, and 
orbital angular momentum states of 4q. States are labeled by 
an index n which differentiates the symmetry type for a given 
parity, the angular momentum J4 of the four quarks, and 
the parity of the state when the strange antiquark is coupled 
TT = —n4, \n; J^). All the states are to be coupled to overall 
complete antisymmetry. States n — 3,4 may be coupled to 
angular momentum J4 = 0,1,2. J4 = 2 is ignored in this 
work. 



to have higher interaction energies than either S = 0,1. 
Color of 4(7 is C = 3 since we require the state of Aqq 
to be a color singlet. The antiquark is an s with spin i, 
spin projection Sz, color c and orbital angular momen- 
tum £,me. The square brackets [ indicate coupling 
of the 4q and q to total angular momentum J, Jz . We will 
consider states of cither parity, tt = ±1 in this work. The 
state in Eq. Q is translationally invariant and requires no 
correction due to center-of-mass motion. 

Consider the fully (pair) correlated state vector j^'s) in 
a five-particle basis of spin, isospin, and color. We have 

2^ 2 3 

l*5(R)) =Y.Y.Y1 i'n^,n,,n,{R)\ns,nt,np) (7) 

ris— 1 nt — 1 rip — l 

where the sums are over the 2^ = 64 spin states Ug, the 
T — isospin space is spanned by two states, nt — 1,2, 
and color states rip = 1,2,3. Tensor interactions in the 
pair potentials v^^'"\rij) populate all of the spin states. 
In the color (isospin) sector we may work in a reduced ba- 
sis since the interactions are mediated by isovector (color 
octet) bosons. The color states of Aq span the C = 3 
irreducible representation (irrep) of SU{3) and the (3, 1) 
irrep of the symmetric group on four objects, 5*4. These 
will be combined with appropriate states in isospin, spin, 
and orbital angular momentum to yield states of com- 
plete antisymmetry. An orthogonal transformation re- 
lates these states to those of Eq.©. 

The OGE and OPE pair potentials, the confining flux 
tube potential, and the two-body correlation operators, 
Fij derived from these all depend, variously, on the state 
\ns,nt,np). Matrix elements of potential pair operators 
appearing in the Hamiltonian in {i9('^) (r,- -) are evaluated 
in Monte Carlo as described in Ref.'29|. The state of 
the gauge field depends on the color configuration of 
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of the pair S = 0, the pair correlation operator satisfies 



TABLE II: Variational energies of the states of Table Q] in 
GeV. The mass is reported as Mg+{n;J2) = (T) + (V") - 
385.5(5) MeV. Statistical errors for Mg+ are < 5 MeV. 



the quarks. Three flux tube topologies are displayed in 
Fig.Q. In the adiabatic approximation, the locations of 
the Y junctions are determined by the minimization of 
the flux tube length. 

The states in Table are obtained by coupling the 
flavor and spin states into irreps of S'4, (q;)tSi where 
T = 0, 5* = 0, 1 and (a) is the irrep, given by the 
Young diagram in the Table. Then couple color and or- 
bital angular momentum sectors {(3)cl, where C = 3 
and L = 0,1. Finally, the flavor-spin {TS) and color- 
orbital (CL) are coupled to give an overall singlet (1*) C 
(a)TS X {(3)cL- The possible combinations of inner prod- 
ucts over spin, isospin, color, and orbital angular mo- 
mentum are given schematically in Table There is 
one negative parity state, |1;1~) and six positive par- 
ity states, |n = 1, • • • ,4; J4). The negative parity state 
|1; 1~) corresponds to the state of 4q with positive par- 
ity, 7r4 = +1 with zero orbital angular momentum. The 
positive parity states, \n = 1, . . . , 4; J^*") correspond to 
states of Aq with negative parity and include one unit of 
orbital angular momentum. We have chosen to neglect 
the state in which the strange antiquark carries one unit 
of orbital angular momentum in the present study. This 
is a prejudice based on diquark correlated wave functions 
[2^. Of course, it only contributes to the higher energy 
positive parity states but a more complete study should 
consider this type of wave function. 



Now that we have the uncorrelated states, we consider 
the correlated state vector, j^E's) of Eq.lO). The determi- 
nation of the ground state of the many-body Schrodinger 
equation can be formulated as a constrained variational 
problem. We wish to minimize (H) subject to the con- 
straint that (^'sl^s) remains constant. This problem 
simplifies greatly working within the ansatz of Eq.jS)) 
but solution of the constrained variational problem now 
yields an upper bound on the ground state energy. The 
ansatz determines two-body equations for the Fij which 
depend on the reduced mass, fiavor, spin, color, and or- 
bital angular momentum of the pair. Pair correlation 
operators for light quark pairs are denoted Fjj'^\ while 

light quark-strange quark pairs are f'/j''''' . For example, 
in the case of two light quarks which have total isospin T, 
color C, orbital angular momentum L, and choosing spin 
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The pair potentials in Eq.(^ are projected into TSC 
channels, wf^^''^ (r)|TS'C) = VTSc{r)\TSC) and the func- 
tions \tsc {r) are generalizations of Euler-Lagrange pa- 
rameters arising in the constrained variational problem. 
Their long distance form ensures that boundary condi- 
tions requ ired of the many-body wave function are sat- 
isfied 30]. They include a small number (^10) of varia- 
tional parameters which are tuned to minimize the vari- 
ational energy. In the case 5 = 1 the tensor interaction 
couples partial waves L to L ± 2 yielding coupled dif- 
ferential equations. The form of the F/j''' is chosen to 
be 
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}(g){]l,T, -Tj}. 



Here Up''-* — fp'^'^^'/fc'^'^^ and the J2p includes hyperfine, 
(Ti ■ (Tj and tensor, Sij operators and isospin. 



and 



color dependent • Tj . The Fj^'"' have only confinement 
and one gluon terms. 

The evaluation of expectation values of the variational 
wave function j^E's) are effected via Monte Carlo integra- 
tion of the 15 dimensional space R = (ri,...,r5). A 
large number (~ 10^) of configurations are sampled from 
the weight function |^'g(R)^'5(R)|. Statistical errors are 
calculated by block averaging over ~ 500 configurations 
[sif . Results for the variational energies of negative and 
positive parity states are shown in Table ITU The mass of 
the state is reported as 

Me+in; J^) = 4m, + + (T> + {V) - Voi^qq) (10) 




FIG. 2: Radial diquark-diquark pair distributions in the state 

|i;i+>- 
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where (T) and (V) are the variational kinetic and poten- 
tial energies, respectively, for each of the states in Table 
in and the states are designated by the symmetry index 
n, the total angular momentum of the four light quarks, 
and the parity tt = —114 of the five-quark state J4 . In 
reality, the true positive parity ground state is some ad- 
mixture of the six positive parity states. However, the 
lowest energy state dominates since the states \n; J4 ) are 
orthonormal: applying the correlations as in Eq.lQ in- 
troduces off-diagonal elements which arc 0{uj?'^\u^^^), 
quite small. States with n 7^ 1 contribute < 10%. 

The constant Vb(4gg) is extrapolated from the single 
baryon spectrum. In the baryon sector the constant 
Vb(3q) is fixed by the nucleon mass niN = 939 MeV 
to be 1312.6(5) MeV where the nucleon mass and sta- 
tistical error were calculated in Monte Carlo with 10^ 
configurations. We assume that Vn scales linearly with 
the number of flux tube ends 18, 32], independent of the 
flavor of the associated quark. This gives for the combi- 
nation Aniq + nris — V{){Aqq) the value —385.5(5) MeV for 
light quark mass vriq = 313 MeV and strange quark mass 
ms ~ 550 MeV. One might take the perspective, as some 
authors have, that the constant should be fitted to the 
mass of the ostensible state 9^ . This would give a large 
value to the constant V()(4gg)/5 = 570 MeV compared 
with the present value Vo(4g^)/5 = 440 MeV. Of course, 
our model neglects effects of chiral symmetry and cannot 
rule out such an effect. We note in this context that the 
present model Hamiltonian does not reproduce the cor- 
rect chiral mass coefficients 33. .34] , a discrepancy that 
is under current study and should be addressed in future 
work. 

Comparison of the CQM with lattice results at large 
pion mass is meaningful since the effects of chiral sym- 
metry are negligible. We note that the negative parity 
state ll;l~) lies below the positive parity state Jill"*"), 
consistent with the lattice data in Refs.j23, |2a]. This is 
true despite the fact that the positive parity state has 
much stronger attraction from OGE and OPE contribu- 
tions, attributable to the high degree of symmetry of the 
TS state (see Table EJ. The {V) of the negative parity 
state is 180 MeV below the lowest positive parity state. 
Exciting one unit of orbital angular momentum on the 
other hand raises the kinetic energy significantly, about 
450 MeV. 

The diquark content of a state gives insight into its 
short-range structure. In Figure El we plot the diquark- 
diquark {DD) pair density as a function of radius, 
for scalar-scalar {SS) p^^^\r) and vector- vector (VV) 

P^DD^i''') diquark configurations for the state ]1; 1+). At- 
traction in the scalar diquark r = 0,5 = 0,C = 3 
SS channel favors this correlation as noted in Ref. [ioj ]. 
It's about 10 times larger than the VV correlation but 
only about 28% of the total strength of the state. Cor- 
relations have a large effect on this non-perturbative fea- 



ture of the dynamics: the uncorrelated state <I'5(1;1"'") 
has equal parts SS and VV, constituting i of the total 
strength. Pair correlations induced by the operators Fij 
in Eq.© therefore significantly increase the strength of 
the SS configurations to those of VV but slightly reduce 
the relative importance of these correlations in the over- 
all wave function. The diquark-triquark correlation of 
Ref . [2lj] is orthogonal to the uncorrelated state ]1; 1+) (it 
appears in states with n 7^ 1) and is strongly suppressed 
in the present calculation. 

We have shown that within the non-relativistic flux 
tube CQM with OGE and OPE interactions the nega- 
tive parity state is lower than the positive parity state 
by ~ 280 MeV when the 0+ is assumed to be a narrow 
resonance. Though the higher symmetry of the lowest ly- 
ing positive parity state significantly decrease its poten- 
tial energy, its unit excitation of orbital angular momen- 
tum raises the energy above that of the negative parity 
state. While diquark structures are a significant factor 
in determining energies of the states, the requisite overall 
antisymmetry of the light quark wave function and pair 
correlations imply that the majority of the wave function 
is described other quark-pair configurations. 
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